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We discuss dynamic spin susceptibility (DSS) in two-dimensional (2D) Dirac electrons with spin- 
orbit interactions to characterize topological insulators. The imaginary part of the DSS appears as 
an absorption rate in response to a transverse ac magnetic field, just as in an electron spin resonance 
experiment for localized spin systems. We found that when the system is in a static magnetic field, 
the topological state can be identified by an anomalous resonant peak of the imaginary part of the 
DSS as a function of the frequency of the transverse magnetic field uj. In the absence of a static 
magnetic field, the imaginary part of the DSS becomes a continuous function of u with a threshold 
frequency lOc- In this case, the topological and the trivial phases can also be distinguished by the 
values of cUc and by the line shapes. Thus the DSS is an experimentally observable physical quantity 
to characterize a topological insulator directly from bulk properties, without observing a topological 
transition. 

PACS numbers: 72.80.Vp, 71.70.Di, 73.43.-f, 76.40.-l-b 

Introduction. Recently, there is growing interest in 
the study of two-dimensional (2D) topological insula¬ 
tors (TIs). This was theoretically predicted by Kane 
and Mele^ based on a model describing electrons on a 
graphene like honeycomb lattice with spin-orbit interac¬ 
tions. Although the first experimental discovery of a TI 
was in HgTe quantum wells,which is described by the 
Bernevig-Hughes-Zhang (BHZ) model,® there are many 
candidates of TIs that have a honeycomb lattice struc¬ 
ture as was originally discussed by Kane and Mele, and 
they are being intensively studied both theoretically and 
experimentally. One of these materials is a silicene,®^^^ 
which is a 2D crystal of silicon. There are also similar 
materials called germanene and stannene, that consist 
of Ge and Sn, respectively.®^’^® There are other types 
of honeycomb lattice materials that consist of two com¬ 
ponents, such as molybdenum dichalcogenides (M 0 S 2 , 

MoSe 2 , etc.).®^’®® These materials have buckled honey¬ 
comb lattice structures with relevant intrinsic spin-orbit 
couplings as compared to graphene. A tunable band gap 
can also be introduced by applying a perpendicular elec¬ 
tric field to the material sheet. 

For these systems, the low-energy electronic properties 
can be described by the 2D Dirac Hamiltonian, as those 
of graphenes, with the Fermi energy at the Dirac point. 

The band gap and the spin-orbit coupling appear as the 
mass term. In such a Dirac system, there is a topologi¬ 
cal phase transition from a TI to a trivial band insulator 
(BI), at a charge neutrality point. The TI has quan¬ 
tized spin Hall conductivity when the spin is conserved, 
and it is more generally characterized by a Z 2 topological 
number.® Since the experiment to observe the topological 
state depends on transport measurements, it is desirable 
to find a non-contact method to identify whether the sys¬ 
tem is a TI or a BI from the bulk properties. 


There are several physical quantities and experiments 
proposed to identify the topological states. For ex¬ 
ample, optical responses,®®“®® spin and valley Hall ef¬ 
fects,dynamical polarization function,^®“^® anoma¬ 
lous spin Nernst effect,^® quantum oscillations, and or¬ 
bital magnetism^'’^®® have been proposed as such ex¬ 
periments. However, these experiments are not simple 
enough to detect TIs, in a sense that they are not direct 
information of the topological state. Especially, in many 
cases, a TI is identified by observing a topological tran¬ 
sition from a BI. Therefore, if the internal parameters 
of the system are not controllable from the outside, the 
identification becomes difficult. 

In this Rapid Communication, we turn our attention to 
the dynamic spin susceptibility (DSS) whose imaginary 
part gives us a very simple index of the topological state 
only by the existence of a certain resonant peak struc¬ 
ture. This method enables us to identify a topological 
phase directly without observing a topological transition. 
We also discuss that an absorption rate in response to a 
transverse ac magnetic field is related to the imaginary 
part of the DSS, just as in an electron spin resonance 
(ESR) measurement. 

2D Dirac fermions. We consider a Kane-Mele type 
Hamiltonian® describing electrons on a honeycomb lat¬ 
tice with an alternating potential A and a spin-orbit in¬ 
teraction K, 

= + X (1) 

(b) * <(b)> 

where c| (ci) is a creation (annihilation) operator at site 
i (spin indices are omitted), {ij) and ((ij)), denote a 
nearest and a next nearest pair, respectively, rji = 1 (rji = 
— 1) for the A (B) sublattice, and I'ij = (2/V3)(di X 
d 2 )z = il, where di and d 2 are unit vectors along the 
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two bonds on which the electron hops from a site j to i. 
Sz is a Pauli matrix describing the electron’s spin. 

In the continuum limit, the local Hamiltonian in a 
static magnetic field is given by 

=hv{^TTxTx + T^yTy) + ATz - K^TzSz, (2) 

where v is the Fermi velocity, ^ = ±1 denotes K and K' 
points, respectively, tt = -VeAjc is the momentum 

operator with V x A = (0,0, H), and A is an energy gap 
related to the alternating potential of the A,B sublattices 
of the honeycomb lattice, (a = x, y, z) are the Pauli 
matrices for the sublattice. 

Because the momentum operators t:± = tt^, ± iiTy fol¬ 
low the commutation relation [7r+,7r_] = — they 

are related to creation and annihilation operators, 
and a, as = \/2ja\ 7r_ = v^fa, for eB > 0 with 
I = \Jch,l\eB\. Then the eigenvalues and eigenstates of 
the Schrbdinger equation for valley-^ and spin-s sector 

n,^^s)) are given using the 

number state \n) of a^a as 




sgn{n)^2fi^v'^\n\/P + (n ^ 0) 

-^A^s (n = 0) 




ll’^l -C-h) 


where A^^ = A — ^sk, = (1 ± ^)/2, and 


(3) 

(4) 


= / ^\J +A^s|/2|i?n I (?T- 7 ^ 0) ^ 

1 (n = 0) ’ 

^ I sgn(n)^lE^'^ - AM^nl (n ^ 0) . 5 ^) 

” 1 e+ (n = 0) ■ 


Here, the Zeeman effect has been ignored, assuming that 
the effective mass of the Dirac system is generally very 
small, but it can easily be taken into account only by 
shifting the Landau levels zL^g/J,BB according to the di¬ 
rections of the spins. The breaking symmetry for positive 
and negative energies at n = 0 eigenvalues is called the 
“parity anomaly” which plays an essential role in the 
topological properties of Dirac systems. Due to the parity 
anomaly, quantized spin Hall conductivity of the system 
at the charge neutrality point with B = 0 becomes 

= ^sgn{K)e{\K\ - |A|). (6) 

This means that the system is TI for |k| > |A| and BI 
for |k| < |A|. 

Properties of DSS. Now let us consider the properties 
of DSS in these two phases of the Dirac system. The 
spin-spin correlation function is given in a Matsubara 
form as 


Xxx{t) 


Xxx(l^m) 


= -(r.5"(r)A"(0))o, (7) 

^Tr G{ek,iujn){ToSx)G{ek,i^n)i'ToSx), 

Ptl 

k,n 

( 8 ) 



FIG. 1: (color online) Imaginary part of the DSS 

—cjim Xxx(u!) (in units of /P) with /i — 0 and F/k = 0.02 
for (a) BI (A/k = 1.2) and (b) TI (A/k = 0.8). The param¬ 
eters are B = 0.1 T, 1 ; = 5.4 x 10® m/s, and k = 4.0 meV 
assnming a silicene. There are regnlar peaks for both cases 
which stem from transitions between —n and n Landan levels 
with |n| > 1. In addition to those, an anomalous peak due to 
n = 0 Landau levels appears at uj = 2^1 h (independent of the 
strength of the magnetic field) only for TI. When the Zeeman 
effect exists, e.g., q^bB/k, = 0.1, the regular peaks split into 
two, while the anomalous peak shifts only to one direction. 


where (• • -/g = Tr • • • e“^^/Tr e“^^, with inverse tem¬ 
perature p. To is the unit matrix, and the temperature 
Green’s function is given as Gispi-On) = — (iw„ — (e — 
/i)/ft -I- isgn(w„)r/Ii)“\ a;+ = -I- Vm with ujn and Vm 

being Matsubara frequencies for fermions and bosons, re¬ 
spectively. The chemical potential and the impurity scat¬ 
tering time are denoted by pL{= 0) and hf2T, respectively. 
After analytic continuation ium —7 w -I- iO, the retarded 
spin-spin correlation function is obtained as 




’ ^,s n,m 




(9) 


where /(e) = -|- 1)“^ is the Fermi distribution 

function, s means the opposite spin of s, and X{x^ y, uj\ H) 
is defined as follows 


■E 


ph ^ (iw„ - xlK)(\Lot - ylh) 


= h 


dfl[— f {hfl)]X{x, y,uj', H). 


( 10 ) 


The transverse spin susceptibility is given by energy tran¬ 
sitions between Landau levels labeled by the same abso¬ 
lute value of n with opposite signs. This selection rule is 
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FIG. 2: (color online) Landau level structures of Dirac 
fermions with alternating potential A and spin-orbit inter¬ 
actions ft for (a) the BI (|ft| < |A|) and (b) the TI (|ft| > |A|) 
states. The two phases are characterized by the configurations 
of the n = 0 Landau levels around the Fermi level (/r = 0). 
For the BI, the n = 0 Landau levels have an opposite sign in 
each spin, so that the quantized Hall conductivity canceled 
out when the magnetic field is tuned off, and the absorption 
rate does not show the anomalous peak, because the transi¬ 
tions between n = 0 Landau levels are not allowed. 


much simpler than that in the current-current correlation 
function for the optical conductivity, which is related to 
the transitions between the |n| and |n± 1| Landau lev- 

els. 17-19.22 

Now we turn our attention only to the imaginary part 
of the DSS multiplied by the frequency —wim Xxxi^^), 
which is experimentally observable as discussed later. 



FIG. 3: (color online) Imaginary part of the DSS, 
—ujlm (in unit of with B — 0 and F/k = 

0.02 for several values of A/k = 0.4-1. 6 . For small F, 
—(ulm Xxx{u)) is discontinuous (continuous) function at the 
threshold frequency ujc = 2ft/ft (ujc = 2A/h) for TI (BI). 


The real part can also be obtained via the Kramers- 
Kronig relation. As shown in Fig. 1, —wim Xxx(j~^) at 
/i = 0 has several peaks corresponding to the transitions 
between the Landau levels below and above the Fermi 
level with the same absolute values of the index |n| > 1. 
We also hnd that an anomalous peak appears only for 
the TI region (|ft| > |A|), whose resonant frequency is 
independent of the strength of the magnetic field. Actu¬ 
ally, for the clean and zero temperature limit (F —^ 0 and 
T —>■ 0), we get the following i5-function peaks, 


Im Xxx{uj) = 


^2(A+l4iJ+i.4-iJ+i.tyl+i.4)2 


■n=l 


X (5(a; - (A+i’^ -b E+^’^)/h) + - A'^)5{u: - 2|ft|/?i) 

( 11 ) 


The last term of Eq. (11) is the anomalous peak at w = 
2k./ h which stems from the transition between the n = 0 
Landau levels with opposite spins. This peak can be 
observed only in the TI where one of the two levels is 
below the Fermi level, as shown in Fig. 2. Physically, 
this peak is related to an energy 2k to make an edge state 
with the spin current in the opposite direction via spin 
flipping. Therefore, we can identify whether the system 
is a TI or a BI, by observing this anomalous peak. The 
numerical value of the frequency of the anomalous peaks 
is estimated as a; ~ 12 THz for silicene with k = 4.0 
meV.^^ 

We can also show that the anomalous peak can easily 
be distinguished from the regular peaks, even if the Zee- 
man effect is taken into account. In the presence of the 
Zeeman effect, Eq. (9) is modified as 
—>■ X{E^^ + sgpL-QB/2,E^ + syfinB/2 ,uj-V.). As shown 
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in Fig. 1, the regular peaks split into two parts Aw = 
while the anomalous peak shifts to one direc¬ 
tion Aw = sgn{K)^BB/h. This is because the valley de¬ 
generacy is lifted by the Zeeman effect for |n| > 1 Lan¬ 
dau levels, while the Landau levels are uniformly shifted 
for n = 0. Therefore, the DSS is still effective informa¬ 
tion to characterize the TI in the presence of the Zeeman 
coupling, if it is sufficiently smaller than the spin-orbit 
coupling, as usually expected for TIs. Furthermore the 
Zeeman shift of the anomalous peak has information to 
identify the sign of the spin-orbit coupling n. 

Case without static magnetic field. So far we have 
assumed the Landau quantization by a static magnetic 
field in a bulk system, so that the peak structures of 
Im Xxx(cj) vanish in the B = 0 cases, and it becomes a 
continuous function of w. However, as shown in Fig. 3, 
the threshold value of the frequency Wc which gives the 
minimum edge of finite —wim Xxxi>x!) is Wc = 2Alh for 
the HI and Wc = 2k/ h for the TI. This difference can be 
easily understood from Eq. (11) for B —>■ 0 which gives 
fiwc = |A+iq-| -I- |A+i j^|. The line shape of Im Xxxi^xj) is 
also different for the two phases: Im Xxx{^) grows contin¬ 
uously from 0 for the BI while it grows discontinuously 
for the TI. This can also be understood by examining 
the i? —>• 0 limit of Eq. (11): For |A| |k|, we get 

-Im Xxx{2.A/h) = 0 and lim^^_^oo[-w Im Xxxieo)] = 
while for |A| <C |k|, we have —Im Xxx{‘2‘H/h) = and 

lim,^^.oo[-w Im Xa:x(w)] = These features may 

also give supplemental information to distinguish a TI 
and a BI. 

Possible experiments. Finally, let us consider how we 
observe Im Xxx(pj) in real experiments. For an example, 
we examine a 2D electron system with a static magnetic 
field and a transverse ac magnetic field. The Hamiltonian 
of the system is expressed as 

'H{t)='H + V{t), (12) 

where % and V (t) are the time-independent and time- 
dependent parts of the Hamiltonian, respectively. For 
the static and the dynamic magnetic fields, 

JB(t) = JB-f J5R,cos(wt), (13) 

with B — (0,0,5) and = (5r, 0, 0), we may choose 
the vector potential as A = 5(—y/2, a:/2,0), Ar, = 
5r( 0, 0,y). In order to discuss general situations, we 
first consider a “non-relativistic” system with a parabolic 
band, and with an electron mass m. Then Eq. (12) is 
given by 


n ^UoiA) + pplbB j drS'^ir), (14) 

V{t) =g/iB5RC0s(wt) J drS^{r) 

- cos (wt) J d^AR(r) ■ Jd(T-), (15) 


where 5“ = / drtplsfgitjjg, with an electron operator ips. 
Jd = ——AB'ipl'ifs means a diamagnetic current induced 
by the dynamic magnetic field, whereas a paramagnetic 
current Jp = — — [V^/)J(r)]^/;s} is included in 

Ho- We have also used the relation ^r(’') • Jp(,r) = 0, 
which is satisfied in 2D systems. 

Then the absorption rate of the dynamic magnetic field 
is obtained up to the second order of V as 

I{u;)=I^^^Tr[p{tW)] (16) 

_ {gpBBR)^u} R 

~- ^ -Imx„^(w), (17) 

where T = 27r/a;, p is a density matrix of Eq. (12), and 
X^x i^) 111® Fourier transform of the retarded transverse 

spin-spin correlation function, 

xi(t) = -i0(t)([5rot(t),5f„,(O)])p. (18) 

Here, we should note that the orbital contributions of the 
ac field do not appear in Eq. (17), because AR(r) • Jd{r) 
is canceled in the time integral, and /(w) becomes just 
the same form as a usual ESR formula for localized spin 
systems.This situation is realized only in 2D systems 
where the electrons do not have a path to move along 
the 2 direction. In Dirac systems, the current operator 
does not have a diamagnetic part because of the linear 
dispersion relation, so that the last term of the right hand 
side of Eq. (15) does not exist. Thus Eq. (17) can also 
be applied to current 2D Dirac systems. 

Summary. We have discussed the dynamic transverse 
spin susceptibility (DSS) of a 2D Dirac system with a 
spin-orbit interaction k and an alternating potential A 
(2). Its imaginary part is related to the absorption rate 
in response to a transverse ac magnetic field. When a 
static magnetic field is applied to the system, the imag¬ 
inary part of the DSS shows an anomalous peak at fre¬ 
quency uj = 2k/ h for the TI. This is related to an energy 
to make an edge state with the spin current in the oppo¬ 
site direction via spin flipping. On the other hand, when 
the static magnetic field is turned off, the imaginary part 
of the DSS becomes a continuous function of lo with dif¬ 
ferent values of threshold frequencies, = 2k/ h for the 
TI and Wc = 2A/fi for the BI, respectively. These prop¬ 
erties enable us to identify the TI directly only from the 
bulk information, without observing a topological transi¬ 
tion from the HI. The DSS is considered to be a relevant 
probe for TIs in more general cases. However, in the 
case of the BHZ model® which describes a TI realized in 
HgTe quantum wells, the present theory cannot be ap¬ 
plied straightforwardly, since the n = 0 Landau levels de¬ 
pend on the strength of the magnetic field^. We also need 
extended discussions for the original Kane-Mele model in¬ 
cluding the Rashba interactions^ where the z-component 
of spins is not conserved. 
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